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ABSTRACT 

We calculate the chiral corrections to Weinberg's prediction for the S-wave iiN scat- 
tering lengths up-to-and-including order 0{M^) making use of heavy baryon chiral 
perturbation theory. For the isospin-odd scattering length a~ these corrections are 
small and bring the prediction closer to the empirical value. In the case of the isospin- 
even a"*" large cancellations appear so that the 0{M^) result depends sensitively on 
certain resonance parameters which enter the calculation of the contact terms present 
at next-to-leading order. 
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1. One of the most spectacular successes of current algebra in the sixties was Wein- 
berg's prediction [1] for the S-wave pion-nucleon scattering lengths, ai/2 = M-,^/4nF^ = 
—2as/2 = 0.175 M~^, with M^^ the physical pion mass and the pion decay constant. 
Tomozawa [2] also derived the sum rule ai/2 ~ <^3/2 = ^^Mt^/SivF^ = 0.263 M~'^. Em- 
pirically, the combination (2ai/2 + 0,3/2)/'^ is best determined. The Karlsruhe-Helsinki 
group gives 0.083 ± 0.004 [3] consistent with the pionic atom measurement [4] of 
0.086 ± 0.004 M~^. The value of ai/2 — 03/2 is more uncertain. The KH analysis leads 
to 0.274 ± 0.005 [5] whereas the VPI group has recently given a larger value [6]. 
Their analysis was critically reexamined by Hohler [7]. In what follows, we will use the 
central values from the work of Koch [3], namely ai/2 = 0.175 and 03/2 = —0.100 
M~^. The agreement of the current algebra predictions with these numbers is rather 
spectacular. However, in the last decade it has become clear that current algebra is only 
the first term in a systematic expansion of the QCD Green functions in powers of (small) 
external momenta and the (light) quark masses [8,9]. Therefore, one would like to know 
what the next-to-leading order corrections to the original predictions are. This is ex- 
actly the question we will address here. The basic framework to perform the calculation 
of these corrections is baryon chiral perturbation theory which makes use of an effec- 
tive Lagrangian of the asymptotically observed fields. In particular, we use the heavy 
fermion formulation [10] in which the nucleons are considered as static sources and one 
has a one-to-one correspondence between the loop and the small momentum expansion. 
We will work out one loop and counter term contributions up-to-and-including order 

based on the power counting scheme developed by Weinberg [8] and later extended 
to the baryon sector by Gasser et al. [11]. For a review on these methods, the reader is 
referred to ref.[12]. 

2. Gonsider the ttN forward scattering amplitude for a nucleon at rest with its four- 
velocity given by = (1, 0, 0, 0).* Denoting by b and a the isospin of the outgoing and 
incoming pion, in order, the scattering amplitude takes the form 

T^" = T+ (a;)5^" + T" (a;)ie^" V= (1) 

with q the pion four-momentum and co = v-q = qq. Under crossing {a b, q ^ —q) the 
functions T+ and T~ are even and odd, respectively, T^{cjij) = ±T'^{—u>). At threshold 
one has q = and the pertinent scattering lengths are defined by 

* Remember that we treat the nucleons as very heavy fields. 
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with m the nucleon mass. The S-wave scattering lengths for the total isospin 1/2 
and 3/2 are related to a"^ via 

ai/2 = a+ + 2a~, 03/2 = a+-a~ (3) 

The abovementioned central empirical values translate into a"*" = —0.83 • 10~^ 
and a+ = 9.17 • 10"^ M'^. In what follows, we will not exhibit the canonical units of 
10-2 M-^. The benchmark values are therefore a+ = -0.83±0.38 and a" = 9.17±0.17* 
compared to the current algebra predictions of a+ = and a~ = 8.76 (using = 138 
MeV and = 93 MeV). 

3. To calculate the scattering lengths, we use the effective pion-nucleon Lagrangian. 
We work in flavor SU(2) and in the isospin limit = nid = rh. The pion fields are 
collected in the matrix U{x) = exp[ir ■ 7r{x)/FT^] ~ u^{x). The effective Lagrangian to 
order 0{q^), where q denotes a genuine small momentum or a quark mass, reads 

r - /-(i) I /-(s) , A3) 



C^^l = H{ivD + gAS-u)H 

~ 8m 



C^^l = ci HHTr{x+) + (c2 - i^)Hv -uvuH + csHu- uH 

^ 8m ^ 

2 

= (^1 - ^) V'v^PHK^xpH + 62 vPHKji^H + 63 v^HKi^^H 



with 



u^x=iv)V^VxUv} (5) 

where H denotes the heavy nucleon field, 5"^ the covariant spin-operator subject to the 
constraint v-S = 0, the covariant derivative acting on the pions and we adhere to the 
notations of ref.[13]. The superscripts (1,2,3) denote the chiral dimension. The lowest 
order effective Lagrangian is of order 0{q). The one loop contribution is suppressed 
with respect to the tree level by thus contributing at 0{q^). In addition, there 
are contact terms of order q^ and q^ with coefficients not fixed by chiral symmetry. 
Notice furthermore that one in addition has to go further in the 1/m expansion of the 
relativistic tree level graphs as can be seen from the terms which come together with 



We have added the errors obtained for ai/2 and 03/2 in quadrature. 
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the ones proportional to ci and h\. At next-to-leading order, all these terms have to 
be retained. Due to crossing symmetry, contributes only to T"'"(a;) whereas C}^1^ 

solely enters T~ {uj). We have not exhibited the standard meson Lagrangian +£itr- 
The contact terms appearing in /^i'^i only contribute to the shift of the pion mass, the 
pseudoscalar coupling and the pion decay constant from their chiral limit to the 
physical values. For the following, we define 

^ = — (6) 



TT 



The calculation of the scattering lengths is straightforward. For the isospin-odd a one 
arrives at 



or = a- (M^) + a- {Ml) + a' (M^) 
a-(M^)=L, a-{Ml) = -Ln 

a- {Ml) = V [1 + f ] + ^ [1 - 2 In ^] - GAttL^M^F^ [b[{X) + h + h] 



(7) 



with A the scale of dimensional regularization. While the constants 62 and 63 are finite, 
bi has to be renormalized as follows to render the isospin-odd scattering amplitude 
T-{uj) finite, 

bi = bl{X)- ^ 



K 



1671^ 



2_p2 

— + -(7.-l-ln4.) 

Here, F is the pion decay constant in the chiral limit. In what follows, we will use 
A = mA = 1.232 GeV, motivated by the resonance saturation principle. Of course, 
physical observables do not depend on this particular choice. Notice that the contact 
term contributions arc suppressed by a factor with respect to the leading current 
algebra term. Matters arc different for the isospin-even scattering length a"*". It consists 
of contributions of order and Ml, 



a+ = a+{Ml) + a+{Ml) 

2r!i - 

8m' 



a+{Ml) = 327:FlL'{c2 + C3 - 2ci - |^) (9) 



a+{Ml) = ^gWM^ - S27TFlL''fx{c2 + cg - 2ci - ^) 

The coefficients ci,2,3 are all finite. From the form of eq.(9) it is obvious that the contact 
terms play a more important role in the determination of a"*" than for a~. 



3 



4. The most difficult task is to pin down the various low-energy constants appearing 
in eqs.(7) and (9). Let us first consider ci,2,3- The coefficient ci can be unambigously 
fixed from the pion-nucleon u-term [13], 



Cl 



4M2 



O-7riv(0) + 



9^ 



647rF2 



= -0.87 ±0.11 GeV" 



(10) 



using qa = 1.26 and cr7riv(0) = 45 ± 8 MeV [14]. To estimate the remaining constants, 
we make use of the principle of resonance saturation [15]. It states that to a high degree 
of accuracy the low-energy constants can be calculated from resonance exchanges by 
integrating out the heavy fields from an effective Lagrangian of the pions chirally coupled 
to the various resonances. In the meson sector, this has been shown to work very well. 
We extend this method to the baryon sector since it is essentially the only method of 
estimating the unknown coefficients. From the meson sector, scalar meson exchange can 
contribute to ci and C3, 

(11) 



Cl - 7;(^3\s = Cl 



Cd 
Cl — 

Crr), 



with 2cd/cm = L^/Ls [15]. The central values for the parameters Cd and Cm given in 
ref.[15] are Cm — 42 MeV and q = 32 MeV, i.e. 2cd/cm = 1.56. However, within 
the uncertainty of L5 and Ls, this ratio can vary between 0.75 and 2.25. In addition, 
intermediate A(1232) states give a contribution to C2 + C3. The general TrAAT-vertex 
can be written as 



= [9,^ -iZ+ ^)7m7.] d'^TT- N + h.c. 



(12) 



where denotes the Rarita-Schwinger field and Z parametrizes the off-shell behaviour 
of the spin-3/2 field.* While Peccei [16] fixed Z = — 1/4, a more recent phenomenological 
analysis of Benmerrouche et al.[17] gives a rather wide range, —0.8 < Z < 0.3. Using the 
empirically well-fulfilled SU(4) coupling constant relation Qt^an = '^9ttn/V^ together 
with the Goldberger-Treiman relation, we can cast the A(1232) contribution to C2 + C3 
into the form 



C2 + Cs L = 



9\ A 



2mA{l + Z) + - Z) 



(13) 



It is mandatory to use here the relativistic formulation of the spin-3/2 field since 



otherwise one would miss a contribution of order to 0+. 
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Clearly, this dependence on Z is one of the major sources of uncertainty in fixing the 
values of the contact terms C2 and C3. Furthermore, there is also a contribution from 
the Roper A'"* (1440) resonance to C2 + C3, 

C2 + C3|^, = -r^T^^^^ (14) 

'^^ 16(m + m*) 

with R— 1 ... 1/4 for Qt^nn* — (1/2 •• ■ 1/4) (/tta^ [18]. In complete analogy, one also has 
a A and N* contribution to the low-energy constants 6^, 62 and 63. At A = mA, it can 
be written as 

6UA)+62 + 63|,_^=-^i 

Finally, other baryon resonances have been neglected since their couplings are either 
very small or poorly (not) known.* It is interesting to note that for Z = —1/4 and 
i? = 1, the iV*(1440) contribution to C2 + C3 and to h\{\) + 62 + ^3 is 4 and 12 per cent 
of the A-contribution, respectively. However, in the latter case the contact terms play 
a much less pronounced role as already discussed. 

5. We now present our numerical results. Consider first the scattering length a~ . Using 
= 138 MeV, = 93 MeV, m = 938.9 MeV, Z = -1/4 and = 1, we find 

a- = (8.76 - 1.29 + 1.69) • 10"^ M"^ = 9.16 • 10"^ M"^ (16) 

where we have explicitely shown the contributions of order M^^, and (and mo- 
mentarily reinstated the units). The total result is in good agreement with the empirical 
value. The largest part of the term comes from the pion loop diagrams, it amounts 
to 1.31 for A = TTT-A- Varying Z within its band of allowed values, a~ varies by ±0.15. 
The dependence on R shows up only in the third digit and is thus irrelevant. The 
one-loop corrections bring the the lowest order value closer to the empirical one. One 
might argue that since the contribution is even larger than the one, the chiral 
series has no chance of converging. However, the important one loop effect can only 
show up at order due to the chiral counting. The two-loop contribution carries an 

* A remark on the p-meson is in order. The chiral power counting enforces a ptttt 
vertex of order q'^ of the form C^pJ-^ = gpT^^T^ {P/j.uiu^ ,u^]) [15]. In forward direction 
the contraction of the p~meson propagator with the corresponding ptttt matrix clement 
vanishes. Therefore, one has no explicit p-meson induced contributions to a~ of order 
and q^. 



+ 



R 

32(m + m*)2 



(15) 
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explicit factor and is therefore expected to be much smaUer. Clearly, one would like 
to perform a calculation beyond 0{q^) as done here, but this is beyond the scope of 
this paper. In the case of the isospin-even scattering length a"*", the situation is much 
less satisfactory. There arc large cancellations between the loop contribution and the 
1/m suppressed kineniatical terms of order and M^. For = 138 MeV, these two 
amount to 0.91 — 0.87 = 0.04 in the conventional units. Therefore, the role of the contact 
terms is even further magnified. In fig.l, we show a"*" as a function of Z for our standard 
input and R = 1. The empirical value of a'^ can be obtained for a small and positive 
value of Z, Z ~ 0.15. As shown in fig. 2, a""" varies also strongly as the ratio 2cd/c^ 
changes. Setting Z = 0, the empirical value results for 2cd/cm = 1.32, not far from its 
central value of 1.56. Clearly, a better understanding of these resonance parameters is 
necessary before one can draw a final conclusion on the accuracy of the chiral expansion 
for a"*". It is, however, gratifying that slight variations of these resonance parameters 
allow one to obtain the empirical value. 

6. To summarize, we have used heavy baryon chiral perturbation theory to calculate the 
corrections up-to-and-including order to Weinberg's lowest order (current algebra) 
predictions for the two S-wave pion-nucleon scattering lengths, a"^. To estimate the 
strength of the various contact terms, we have made use of the principle of resonance 
saturation which is known to work very accurately in the meson sector. The main results 
of this investigation are: 

• The chiral corrections to the isospin-odd scattering length a~ are small and positive 
and move the lowest order prediction closer to the empirical value. The main effect 
comes form the pion loop diagrams. The contact term contribution is relatively 
small, thus masking the uncertainty in estimating the coefficients which appear 
together with these terms. 

• The situation is very different for the isospin-even scattering length o"*" . The contact 
term contribution completely dominates the chiral expansion since there is a large 
cancellation between the one-loop and the kinematical corrections. The total result 
for a"*" is very sensitive to some of the resonance parameters, the empirical value of 

can, however, be obtained by reasonable choices of these. 

Evidently, further investigations have to go into two directions. First, a calculation 
beyond 0{q^) has to be performed to find out how fast the chiral expansion of the 
scattering lengths converges. Second, a better understanding of the coefficients of the 
contact terms appearing at order and higher is necessary to further pin down the 
prediction for a"*". We hope to come back to these topics in the near future. 
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Figure Captions 

Fig.l The scattering length as a function of Z in units of 10~^M~-^. The input is 
specified in the text. The empirical range is also shown (the solid line gives the 
central value). 

Fig.2 The scattering length as a function of 2cd/cm in units of 10~^ The input 

is specified in the text. For notations, see fig.l. 
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